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NOVEL PHARMACOPHORE FOR THE DISCOVERY AND TESTING OF NA, 
K-ATPASE INHIBITOR COMPOSITIONS AND METHODS FOR THEIR USE 
IN TREATING CARDIOVASCULAR DISEASES AND CONDITIONS 



CROSS REFERENCE TO RELATED APPLICATIONS 
[0001] This application is a continuation-in-part of United States Provisional Application 
Serial Number 60/425,037, filed November 7, 2002. 

FIELD OF THE INVENTION 
[0002] The present invention relates to the discovery of a three dimensional 
pharmacophore based on a model of the Na, K-ATPase that will enable the discovery and 
design of novel positive inotropic agents using molecular modeling, pharmacophores, rational 
drug design techniques, and methods of using the novel compounds in treating heart diseases. 

BACKGROUND OF THE INVENTION 
[0003] Each of the references cited herein is incorporated by reference in its entirety. 
[0004] Heart diseases, such as cardiac failure, paroxysmal atrial tachycardia, and 
caridomyopathy, are a major cause of morbidity and mortality around the world. Cardiac 
failure occurs when the heart becomes unable to pump blood effectively at a rate that meets 
the needs of the metabolizing tissues. The underlying cause of cardiac failure lies in the 
reduced contractility of heart muscles. Reduced contraction of the heart leads to reduced 
heart output of blood volume while the amount of blood being returned to the heart remains 
the same. The result is an increase in heart blood volume, which in turn leads to heart 
congestion. Congestive heart failure causes lowered blood pressure, poor renal blood flow 
with concomitant edema in the lower extremities and pulmonary edema, as well as renal 
failure. Positive inotropic activity, meaning increasing the force of the contraction of the 
heairt, is thus a primary goal in the treatment of heart failure patients. 
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[0005] While recent advances in diagnostic and surgical procedures have revolutionized 
cardiovascular medicine, similar progress has not been attained with respect to 
pharmacological agents. Increasing the force of contraction of the failing heart is still largely 
achieved by cardioglycosides (also called cardioactive glycosides) derived from the plant 
Digitalis purpurea, (foxglove) such as digoxin and digitoxin. Unfortunately, digoxin, 
digitoxin, and other cardioglycoside drugs, the most widely used clinical inotropic agents, 
have narrow therapeutic indices and a high potential for considerable toxic effects including 
ventricular arrhythmias. Therefore, the difference between the dosage leading to a therapeutic 
increase in heart muscle contraction and the dosage leading to fatal toxicity is minute. This 
potentially deadly property severely curtails the clinical usefulness of the currently available 
drugs. While much effort has been directed toward the discovery of novel therapeutic agents, 
the lack of three-dimensional structural coordinates for the receptor of inotropic drugs, 
Sodium Potassium (Na, K) ATPase, has markedly limited success. Increasing the therapeutic 
index of positive inotropic Na, K-ATPase (also referred to as Na+, K+-ATPase herein) 
inhibitors would thus represent a major advance in the clinical management of heart disease. 
[0006] Sodium Potassium ATPase is a heterodimeric transmembrane protein that actively 
exchanges sodium and potassium ions across cell membranes using energy from the 
hydrolysis of ATP. Inhibition of the Na, K-ATPase in the myocardial cells of the heart 
muscle results in the sequential increase of cytosolic sodium and calcium concentrations, 
thereby increasing the force of contraction. Regulation of the levels of sodium and potassium 
contribute to many essential physiological processes including maintenance of the membrane 
potential for muscle contraction. 

[0007] Cardioglycoside inhibitors bind to the extracellular domain of the catalytic la 
subunit of the Na, K-ATPase and inhibition of the ATPase results in both the beneficial and 
detrimental effects associated with cardioglycoside administration. While researchers have 
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obtained low resolution crystal structures (~6A) of the Na, K-ATPase and gained structural 
insights from solid-state NMR studies of ATPase-inhibitor binding, the dearth of detailed 
structural information has hindered the discovery of novel inotropic inhibitors. One of the 
main obstacles to Na, K-ATPase crystallization is that the protein is structurally complex, 
possessing ten transmembrane spanning domains (H1-H10), which makes the development of 
the high quality crystals of integral membrane proteins required for X-ray crystallography 
exceedingly difficult. 

[0008] The exact mechanism of inhibition has remained elusive. Although researchers 
have examined alternative inhibitory agents, of the hundreds of cardioglycosides tested, none 
have been shown to have larger therapeutic indices than existing therapeutic agents do. 
Moreover the lack of structural knowledge of the Na, K-ATPase has thwarted the structural 
investigation and structure based discovery of novel alternative treatments. A novel approach 
to the design, discovery, and use of inotropic drugs with reduced toxicity would thus be very 
usefiil in the treatment of heart failure. 

SUMMARY OF THE INVENTION 
[0009] The present invention creates a novel pharmacophore from which to generate 
scaffolds, and then molecules, which have inotropic properties. Preferably, the 
pharmacophore is defined by the parameters of Table 4 and Table 5, which reference the 
spheres in Figure 4. The invention also contemplates the method of creating the novel 
pharmacophore. 

[0010] Another aspect of the present invention sets forth a basic scaffold created using 
the novel pharmacophore. The basic scaffold molecule is a novel Na, K-ATPase ligand. The 
scaffold can be used to develop positive inotropic drugs that inhibit Na, K-ATPase activity. 



3 



WO 2004/043384 PCT/US2003/035636 
[001 1] A third aspect of the present invention discloses further developed potential 
scaffolds of the novel Na, K-ATPase ligands. These scaffolds represent the parent molecules, 
which satisfy the basic three point pharmacophores. The Na, K-ATPase inhibitor compounds 
are of the formula: 



O 




wherein Rl, R2, R3 and R4 can be any organic functional group containing a hydrogen bond 
donor or a hydrogen bond acceptor and X is any element or group that allows the compound 
to retain inotropic activity. Preferably, X is N, O, S, or C. 

[0012] A further aspect of the present invention takes the scaffold molecules and 
develops them into functional drugs for Na, K-ATPase inhibition, which can be used to 
prevent or treat heart conditions, particularly congestive heart failure and supraventricular 
arrhythmia. Eight compounds have shown promise as inotropic agents, with three of the 
eight demonstrating high therapeutic potential based on their ATPase inhibitory 
characteristics. 

[0013] These and other aspects of the invention will be further elucidated in the following 
descriptions and examples. The descriptions in the present invention are provided only as 
examples and should not be understood to be limiting on the claims. Based on the 
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description, a person of ordinary skill in the art may make modifications and changes to the 
preferred embodiments, which does not depart from the scope of the present invention. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0014] Figure 1. This figure shows how the threaded alignment of the loc-subunit of 
SERCAla (skeletal muscle sarcoplasmic reticulum/endoplasmic reticulum Ca2 + -ATPase) 
[SEQ ID NO: 1] with the sheep Na, K-ATPase [SEQ ID NO: 2] was accomplished. The 
alignment shown is the result of "threading" the Na + , K + -ATPase sequence to that of the 
template. Identical residues between the two sequences are illustrated with filled black 
boxes, while similar residues correspond to the filled gray boxes. Regions modeled are noted 
by lines above lettered designations, with the letter shown at the beginning and end of the 
sequence encoding the region. The letters correspond to letters attached to the helices and 
loops depicted in Figures 2A and 2B. The solid lines under the sequence of Figure 1 represent 
helical structures, while the broken lines are indication of loop regions. Local areas of high 
identity correlating with functional motifs are shown by open black boxes and described in 
the text. Positions of sequence changes between the sheep and rat Na + , K + -ATPases are 
denoted with an asterisk (*). 

[0015] Figure 2. This figure shows the resulting homology model of the extracellular 
loops and the transmembrane domains regions of the sheep lcc-subunit of the Na, K-ATPase. 
Shown is a visual representation of the three-dimensional coordinates assigned to the 
homology model (2A), and a representation of the cc-helices from the extracellular domain 
(2B). Only alpha caibons are depicted. The numbering, showing the orientation of the 
helices, is consistent with the information obtained from low-resolution structures. 
Tryptophan and lysine residues are displayed (2A). The lettering corresponds to the amino 
acid sequence portions shown in Figure 1. 
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[0016] Figure 3. This figure shows the binding orientation of cardioglycosides to the la 
subunit of the Na, K-ATPase predicted by the methods of the invention. White numbers in 
the figure indicate transmembrane domains. The extracellular (3 A) and cross-section (3B) 
represent the consensus orientation for all nine docked cardioglycosides. Both Na* K + - 
ATPase and ligands are depicted with Connelly surfaces. The figure depicts consensus 
orientation of 25 individual dockings of ouabain to the la-subunit of the sheep Na + , K + - 
ATPase (3C). Orientation of digoxin (3D, upper) and digoxigenin bisdigitoxose (3D, lower) 
illustrate the identical conformation regardless of the number of sugar moieties. 
[0017] Figure 4. This figure shows the pharmacophore elements fitting the binding 
criteria for la-subunit Na, K-ATPase inhibition. Each pharmacophore element as described 
herein is represented by a point encased in a sphere. The spheres (diameter = 2.5A) represent 
the tolerance of location allowed during database searching. The pharmacophore elements 
are combined to form three, four and five point pharmacophores representing increasing 
stringency of search criteria (Table 4). 

DETAILED DESCRIPTION OF THE INVENTION 

Definitions 

[0018] The present invention can best be understood in light of the following 

definitions. 

[0019] "Inotropic" means affecting the force or energy of muscular contractions. 

[0020] "Glycosides" are any compound that contains a constituent sugar in which the 
hydroxyl group attached to the first carbon is substituted with an alcoholic, phenolic, or other 
group. 
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[0021] "Rational drug design" means modeling the molecular structure of the target of a 
drug and then designing a drug that will attack the target 

[0022] "Pharmacophore" is a model for developing one or more molecular scaffolds or 
structure used as the basis for drug development 

[0023] "Scaffold" means a supporting framework, such as the basic backbone or structure 
of a molecule with designated positions capable of containing a one or more molecules or 
functional groups. 

[0024] "Dynamic programming" is a use of scoring matrices and gap penalties to produce 
optimal alignment. 

[0025] 'Unaltered" means not changed, substituted, deleted, or mutated. 
[0026] "Pharmaceutically effective dose" as used herein refers to the amount of an novel 
inotropic compound or novel compound composition comprising one or more of the novel 
compounds described herein that produces a desired therapeutic effect, such as treating the 
target disease. The precise amount of the pharmaceutically effective dose of a novel 
compound or novel compound composition that will yield the most effective results in terms 
of efficacy of treatment in a given subject will depend upon the activity, pharmacokinetics, 
pharmacodynamics, and bioavailability of a particular inotropic compound, physiological 
condition of the subject, including age, gender, disease type and stage, general physical 
condition, responsiveness to a given dosage and type of medication, the nature of 
pharmaceutically acceptable carrier in a formulation, and a route of administration, among 
other potential factors. Those skilled in the clinical and pharmacological arts will be able to 
determine these factors through routine experimentation consisting of monitoring the subject 
and adjusting the dosage. Remington: The Science and Practice of Pharmacy (Gennaro ed. 
20 th edition, Williams & Wilkins PA, USA) (2000). 
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[0027] As used herein, the "target disease" may be any disease in which inhibiting Na, K- 
ATPase treats the disease. Without limitation, the target diseases are heart diseases, 
particularly congestive heart failure and supraventricular arrhythmia. While it is possible for 
a novel compound to be administered as a pure or substantially pure compound, it is 
preferable that the novel inotropic compound be administered as a composition in the form of 
pharmaceutical formulations or preparations suitable for a particular administration route. A 
novel compound composition comprises one or more novel inotropic compounds and a 
pharmaceutically acceptable carrier. 

[0028] The term "pharmaceutically acceptable carrier" as used herein means a 
pharmaceutically-acceptable material, composition or vehicle, such as a liquid or solid filler, 
diluent, excipient, solvent or encapsulating material, involved in carrying or transporting one 
or more novel inotropic compounds from one tissue, organ, or portion of the body, to another 
tissue, organ, or portion of the body. The novel inotropic compounds may be administered in 
liposomes. Each component must be "pharmaceutically acceptable" by being compatible 
with the other ingredients of the formulation. It must also be suitable for use in contact with 
the tissue or organ of humans and animals, specifically the heart and circulatory system, 
without excessive toxicity, irritation, allergic response, immunogenecity, or other problems or 
complications, commensurate with a reasonable benefit/risk ratio. 
[0029] A "route of administration" for a novel compound or composition can be by any 
pathway known in the art, including without limitation, oral, enteral, nasal, topical, rectal, 
vaginal, aerosol, transmucosal, transdermal, ophthalmic, pulmonary, and/or parenteral 
administration. A parenteral administration refers to an administration route that typically 
relates to injection. Parenteral administration includes, but is not limited to, intravenous, 
intramuscular, intraarterial, intrathecal, intracapsular, infraorbital, intracardiac, intradermal, 
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intraperitoneal, transtracheal, subcutaneous, subcuticular, intraarticular, subcapsular, 

subarachnoid, intraspinal, via intrasternal injection, and/or via infusion. 

[0030] 'Treatment" of or "treating" a disease may mean preventing the disease by causing 

clinical symptoms not to develop, inhibiting the disease by stopping or reducing the 

symptoms, the development of the disease, and/or slowing the rate of development of the 

disease, relieving the disease by causing a complete or partial regression of the disease, 

reducing the risk of developing the disease, or a combination thereof. 

[0031] The tenn "contacted" when applied to a cell, tissue or organ means the process by 

which a novel compound or compound composition is delivered to the target cell, tissue or 

organ, or placed in direct proximity of the cell, tissue, or organ. 

[0032] "Therapeutically effective amount" is the amount of novel compound or 

composition that, when administered to a subject, is effective to bring about a desired effect. 

In this case, that effect is typically an antagonistic effect that ultimately decreases the activity 

oftheNa,K-ATPase. 

[0033] 'Thaimaceutically acceptable carriers" for the novel inotropic compositions may 
include sugars, starches, cellulose, powdered tragacanth; malt; gelatin; talc; excipients, oils, 
glycols, esters, agar; buffering agents, and other non-toxic compatible substances employed 
in pharmaceutical formulations. The formulations may conveniently be presented in unit 
dosage form and may be prepared by any methods well known in the pharmaceutical arts. 
Formulations suitable for oral administration may be in the form of capsules, cachets, pills, 
tablets, lozenges, powders, as granules, bolus, electuary, or a paste, as a solution or a 
suspension in an aqueous or non- aqueous liquid, as an oil-in-water or water-in-oil liquid 
emulsion, as an elixir or syrup, or as pastilles, each containing a predetermined amount of one 
or more of the novel inotropic compositions as an active ingredient. 



9 



WO 2004/043384 PCTYUS2003/035636 
[0034] In solid forms for oral administration, capsules, tablets, pills, powders, granules and 
the like, may by used. In the case of capsules, tablets and pills, the pharmaceutical 
compositions may also comprise buffering agents. Compressed tablets may be prepared using 
binder, lubricant, inert diluent, preservative, or disintegrant. Liquid dosage forms for oral 
administration include phannaceutically acceptable emulsions, microemulsions, solutions, 
suspensions, syrups and elixirs. Formulations for the topical or transdermal administration of 
a novel inotropic composition include powders, sprays, ointments, pastes, creams, lotions, 
gels, solutions, patches and inhalants, such as by aerosol. 

[0035] Formulations suitable for parenteral administration comprise a novel inotropic 
composition in combination with one or more phannaceutically acceptable sterile isotonic 
aqueous or nonaqueous solutions, dispersions, suspensions or emulsions, or sterile powders 
which may be reconstituted into sterile injectable solutions or dispersions. The formulations 
may contain antioxidants, buffers, bacterostats, solutes that render the formulation isotonic 
with the blood of the intended recipient or suspending or thickening agents. Formulations 
suitable for parenteral administration may also contain adjuvants such as preservatives, 
wetting agents, emulsifying agents and dispersing agents, various antibacterial and antifungal 
agents, or agents that delay absorption such as aluminum monostearate and gelatin. 

Aspects 

[0036] A novel pharmacophore is disclosed. The pharmacophore allows the generation 
of basic scaffold molecules, which can then be transformed into more detailed scaffold 
molecules having inotropic properties. The pharmacophore is created to generate Na, K- 
ATPase ligands and is defined by the parameters of Table 4 and Table 5 and the spheres in 
Figure 4. The more detailed scaffolds can then be used to create potential drugs, which can 
also be tested within the pharmacophore model for binding potential. The scaffolds at least 



10 



WO 2004/043384 



PCTAJS2003/035636 



satisfy the basic three point pharmacophores, but may satisfy the four or five point 
pharmacophores. 

[0037] The Na, K-ATPase inhibitor scaffold compounds may be of the formula: 



for example, wherein Rl, R2, R3 and R4 can be any organic functional group containing a 
hydrogen bond donor or a hydrogen bond acceptor and X is any element or group that allows 
the compound to retain inotropic activity. Preferably, X is a single element, such as N, O, S, 
or C. 

[0038] The scaffold molecules are the basis functional drugs for Na, K-ATPase 
inhibition, which can be used to prevent or treat heart conditions, particularly congestive 
heart failure and supraventricular arrhythmia. Eight compounds have been tested shown as 
inotropic agents, with three of the eight drugs have demonstrated particularly high potential. 
These drugs are shown below: 

(1) (+>Catechin Hydrate (experimentally noted as Compound 9) 




o 




o 



HO 




OH 
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(2) 4 f -5' Di nitroflorescein (experimentally noted as Compound 14). This compound 
demonstrates strong ATPase inhibition. 




4\5' Dinitroflorescein 

(3) Emodin (experimentally noted as Compound 15). 



OH 


0 


OH 
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(4) Lacmoid (experimentally noted as Compound 19). This compound demonstrates strong 
ATPase inhibition. 




#19 R= NH a or OH 



(5) R-S Dobutamine (experimentally noted as Compound 25). 
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(6) Compound 26 




R1 »OHorNH 2 

(7) Compound 27 




(8) Compound 28, This compound demonstrates strong ATPase inhibition. 




The data for the tests of the eight compounds is found in Tables la and lb. 
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Table la 
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Table lb 



Activity of NKA 


assay started by addition of KCI 


Cardofllycosides will show qreater affinity with KCI added last rather than ATP 
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[0039] The drugs found to be efficacious in treating heart diseases will be administered to 
individuals in need thereof in phaimaceutically effective amounts for the duration of the 
disease, until the disease has abated, or as needed to prevent the disease. The novel inotropic 
compounds should also be administered in a pharmaceutically acceptable carrier. There are 
many methods and routes of administration and possible carriers, which are discussed in 
detail in the definitions section. An ordinarily skilled clinician will be able to determine the 
appropriate amount of the one or more novel inotropic compounds to administer to an 
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individual. Preferably, the novel inotropic compound has a wider therapeutic index that the 
commonly used cardioglycoside drugs, such as dioxin or ouabain, 
[0040] Finally, a general method of making and using the pharmacophore to create an 
Na, K-ATPase inhibitory compound comprises: 

(a) creating alignment between SERCAla and Na, K-ATPase, wherein SERCAla acts as the 
structural template, and the alignment comprises use of both dynamic programming and 
threading techniques; 

(b) transferring coordinates from a template to a model for structurally conserved regions; 

(c) generating variable regions; 

(d) refining the model through energy minimization steps; and 

(e) performing docking analysis of prospective drug candidates. 
The process may further comprise the steps of: 

(f) delineating the essential pharmacophore elements for high binding affinity; 

(g) searching databases of known compounds using the restraints as implicated by the 
pharmacophore with allowable tolerances (e.g. ±10%); and 

(h) utilizing de novo rational drug design and computer aided molecular modeling to design 
novel compounds using the restraints as implicated by the pharmacophore with allowable 
tolerances (e.g. ±10%). 

[0041] As now explained in further detail, the present invention has found a way to 
utilize the recent solution of the structural coordinates of the E* conformation of SERCAla 
(skeletal muscle sarcoplasmic reticulum/endoplasmic reticulum Ca 2 + -ATPase) as a structural 
template for homology three-dimensional modeling of the Na, K-ATPase. SERCA also 
occurs in the E2 conformation, and in SERCAlb, and SERCA2 forms. The Na, K-ATPase 
and other P 2 -type ATPases are thought to share important structural features with SERCAla 
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and sequence analysis of the family members established that these structures are 
evolutionarily related 

[0042] Accordingly, by the methods of the instant invention, a model of the extracellular 
and transmembrane domains of the human la-subunit of the Na, K-ATPase has been 
developed, starting from the high-resolution crystal structure of SERCAla. The first and 
most important step in any modeling process is the alignment of the protein sequence to be 
modeled to that of the structural template. With a sequence identity of only -25%, aligning 
the la-subunit of the Na, K-ATPase to SERCAla by virtue of residue identity/homology is 
thwarted with complications and thus such an afignment cannot be confidently used as a basis 
for a three-dimensional structure. 

[0043] Therefore, in one preferred embodiment, the invention takes advantage of the 
known structural and functional similarities of the family of P 2 -Type ATPases and uses a 
technique that essentially "threads" the sequence of interest to libraries of known protein 
folds from resolved structures, m this manner, the possible tertiary structure configurations 
that provide environments that satisfy the primary sequence of the la-subunits of the Na, K- 
ATPase can be investigated. Through the independent alignment of the sheep la -subunit 
sequences to SERCAla, the instant invention has achieved alignments that have a high level 
of confidence. Areas of high local sequence similarity (indicated by shaded black and gray 
boxes, Figure 1) are well aligned as exemplified by the perfect alignment of the 7 residues 
(DKTTGTLT) (residues 351-357 ofSEQ ID NO:l) that comprise the E1-E2 ATPases 
(aspartate) phosphorylation site that is ubiquitous to P 2 -type ATPases. The TGES (residues 
181-184 of SEQ ID NO:l) sequence prior to the phosphorylation site and the GDGXNDXP 
(where X is any residue) (residues 701-708 of SEQ ID NO:l) motif (Figure 1, open black 
boxes) are further evidence of the alignment. These data demonstrate that after aligning the 
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sequences without bias to residue identity, the resulting orientation effectively aligns these 
essential regions. 

[0044] Further extrapolating from the structural conservation of P 2 -type ATPases, the 
present invention uses computational homology modeling techniques to derive a three- 
dimensional structure for the extracellular and transmembrane domains of the la subunits of 
both human and sheep Na, K-ATPases. Thus, in another preferred embodiment of the 
invention, these three-dimensional models and specific mutations within the extracellular 
domain as they relate to ouabain affinity/binding are applied to examine the mode of 
inhibitor-receptor association. As ouabain and the other cardioglycosides bind to the 
extracellular domain of the la-subunit, five structurally conserved regions (SCRs) gre 
modeled which collectively consist of the ten transmembrane domains and connecting 
extracellular loop regions. These SCRs are notated by solid numbered lines in Figure 1 and 
the numbers correlate with the helices depicted in Figure 2. Areas with extended regions of 
hydrophobicity in the primary sequence are then compared to the transmembranes predicted 
by the alignment, those predicted by secondary structure predictors, and the transmembrane 
domains proposed as a result of hydrophobicity studies. The results of the analyses are 
tabulated (Table 2) and the high level of consistency provide additional evidence to 
corroborate the alignment scheme of the invention. Furthermore, when a standard pairwise 
alignment is performed using a hydrophobicity matrix, the fundamental alignment (Figure 1) 
is maintained. Consequently, supportive data from three independent methods as well as 
structural motif conservation are obtained, all of which converge on the final sequence 
alignment of the present invention. 
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Table 2. 
Helices 



Alignment 



Hydrophobicity Plots 



Meinsat 
89-110 
123-142 
287-305 
314-338 
771-791 
800-817 
842-866 
914-930 
946-964 
978-994 



1 

2 
3 
4 
5 
6 
7 
8 
9 

10 



83-109 
126-149 
285-298 
316-337 
768-789 
799-819 
845-871 
909-934 
948-969 
981-997 



89-110 
123-140 
282-304 
311-342 
770-791 
796-818 
841-865 
909-930 
945-970 
975-996 



[0045] In another preferred embodiment of the invention, the structural coordinates of the 
helical regions of sheep la-subunit are assigned, and loops are generated that are sterically 
acceptable and that result in reasonable phi and psi angles for the spliced regions (Ihsightn, 
Accelrys). The resulting structure is then subjected to a series of refinement steps. Initially, 
the entire protein is relaxed and thefmodynamically energy niinimized to reduce initial steric 
bumps resulting from residue replacement within the secondary structural features. Next, the 
residue side chains are checked for appropriate rotomer positioning using a standard rotomer 
library defined from known crystal structures. Finally, the model is energy niinimized 
through a series of steps in which hydrogens, side chains and finally, the entire molecule are 
successively allowed freedom of movement. Ribbon representations of the cross-section and 
the extracellular side views of the final sheep model are shown in Figures 2A and 2B 
respectively. 

[0046] The resulting model can be structurally validated by the use of the wealth of 
biological data available for the sheep la-subunit for structural validation. Site directed 
chemical labeling of the extracellular loops of the sheep Na, K-ATPase has been utilized to 
suggest structural orientation. Cysteine residues for PI 18, T309, L793, L876, and M973 
have previously been introduced to investigate membrane topology. All of the introduced 
cysteine residues could be chemically labeled indicating that these residues are exposed to the 
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extracellular surface and our model places all five residues in extracellular loop regions. 
Studies have shown that the preference of tryptophan residues and to a lesser extent lysine 
residues for membrane surfaces is a feature ubiquitous to membrane proteins. Analysis of the 
positioning of these residues in the model of the present invention shows that in fact these 
amino acids do tend to associate close to the proposed transmembrane surfaces, providing 
further validation of the model. Moreover, as expected, there appears to be a propensity for 
tryptophan residues to be localized on the extracellular surface while lysine residues appear 
to be incorporated on the intracellular membrane. In summary, both the alignment and model 
of the present invention are strongly supported by experimental data. 
[0047] In another preferred embodiment, with the information obtained from the three 
dimensional structures of the transporter, the invention provides a pharmacophore model for 
the Na, K-ATPase, which describes the types of atoms or pharmacophoric elements and their 
geometric arrangements that are common across the set of active cardioglycosides in relation 
to the constraints imposed by the molecular binding site. The three-dimensional model of the 
pharmacophore of invention makes it possible for the first time to examine the molecular 
mechanism of Na, K-ATPase inhibition, to screen for novel inotropic inhibitors of the Na, K- 
ATPase with higher therapeutic indices, as well as to design such inhibitors through the use 
of rational drug design techniques. With the 3D-model of the Na, K-ATPase provided by the 
present invention, the mode of cardioglycoside binding can now be elucidated. Three- 
dimensional quantitative structural-activity relationship models (QSAR) based on the 
inhibition of the sheep la-subunit as a result of incubation with cardioactive compounds 
previously developed and the resulting comparative molecular field analysis (CoMFA) 
models provide structural implications for the binding site based on the correlation of 
changes in ligand conformation and biological activity. In combination with the model of the 
present invention, these models suggest that the length of the binding pocket is -20 A, which 
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is corroborated by earlier estimates of ~19A. Mutational studies suggest that the ligand 
binding site is situated between H1-H2 (Ribbon No. 1, Figure 2 A) and H5-H6 (Ribbon No. 3, 
Figure 2A), and H1-H2 and H9-H10 (Ribbon No. 5, Figure 2A) corresponding to distances of 
20A and 22A respectively. Collectively these data suggest that the binding site consists of a 
~20A groove comprised of the extracellular loops connecting H1-H2, H5-H6, and H9-H10. 
[0048] Thus, in another preferred embodiment, the present invention provides a three- 
dimensional model of the Na, K-ATPase cardiotonic binding site. In a further preferred 
embodiment of the invention, a pharmacophore model is provided, which precisely defines 
the mechanism of interaction between the Na, K-ATPase and an inhibitory molecule. This 
model allows for the de novo design of drug candidates based on the Na, K-ATPase receptor 
structure. The pharmacophore allows for the development of inotropic agents that are 
structurally unique from those chemical classes known and examined to date, thus making it 
possible to develop inotropic agents with higher (also called wider) therapeutic indices. 
[0049] The systematic docking of 16 cardioglycosides inside the putative Na, K-ATPase 
binding pocket provided by the instant invention was achieved by using ) GOLD (Genetic 
Optimization for Ligand Docking) Jones et al JMolBiol 254, 43-53 (1995). However, any 
accurate docking program is contemplated, including conformational sampling-based 
methods (e.g., free energy perturbation, linear interaction energy approximations), 
knowledge-based "potential of mean force" methods (Pm£ Drugscore), force-field methods 
(GOLD, Dock, AutoDock), Apropros, CombiBtULD, Escher, FlexiDock, FLOG, GRAMM, 
HEX, HotDock, LIGIN, PUZZLE, STALK, and empirical free-energy scoring functions 
(Ludi, Chemscore, Score, Fresno, FlexX, Pip). In addition, the docking and scoring methods 
of US Provisional Application Number 60/477,714, entitled "Novel Knowledge-based 
Docking and Scoring Algorithm for Virtual Drug Screening (Welsh, et al.), filed June 10, 
2003 is incorporated by reference herein. 



20 



WO 2004/043384 PCTYUS2003/035636 
[0050] GOLD is a ligand-docking program, which uses a genetic algorithm (GA) to 
explore ligand conformation and satisfy ligand-binding requirements 
(http://www.ccdc.cam.ac.uk/prods/gold). One advantage over other docking systems is that 
GOLD allows not only for ligand flexibility, but also for limited flexibility of the binding 
pocket creating a more realistic environment for inhibitor association. Further, the stochastic 
nature of GOLD ensures that the search space is well explored and local extremes are less 
influential as compared to deterministic algorithms. Thus, it was possible to select a subset 
of inhibitors consisting of all molecules 15-30A in length and including ouabain, digoxin and 
digitoxin. These molecules best fit the binding site of the invention, and as a result are 
unlikely to have more than one mode of association. The combination of experimental and 
computational results suggests that the ligand binding domain is ~20A which effectively 
defines a docking sphere of a ~10A radius. A centroid of the sphere was determined as the 
intersection of the vectors defining H1-H2 with H5-H6, H1-H2 with H9-H10 and H5-H6 with 
H9-H10. Each inhibitor was docked 25 independent times. 

[0051] To explore more thoroughly the protein-inhibitor interaction, the radii of the 
docking spheres were varied (with increasing radii 10A, 12A, and 15A). Similar results were 
obtained with the two smaller docking spheres while the excessively large volume of the 15 A 
docking sphere resulted in unrealistic interactions in which the molecules docked 
perpendicular to the molecular surface of the Na, K-ATPase. Using the results generated 
with the 12A-docking sphere, the consensus mode of binding for each inhibitor from the 
family of 25 dockings was determined. As all the binding modes were rank ordered using the 
GOLD fitness function, which measures steric and electrostatic complementarity between 

each ligand conformation and the receptor, the molecule corresponding to the highest scoring 

i 

conformation of each binding consensus was selected. As an example, the alignment of the 
25 conformations of ouabain (Figure 3 A) and the overall alignment of all 16 inhibitors 
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representing the consensus conformation (Figure 3B) is provided. It is notable that while 
GOLD is a stochastic algorithm, the independent dockings converged to a single consensus 
both within dockings of an independent compound as well as across dockings between 
compounds. In a manner consistent with experimental findings, the cardiotonic steroid 
moiety and the unsaturated lactone ring are oriented toward H1-H2, while the carbohydrates 
are directed toward H9-H10. Compounds with more than one sugar see these moieties wrap 
around H7-H8. All compounds completely occupy the potential binding pocket which is 
situated over the channel (depicted by a fast Connolly channel surface) thus effectively 
inhibiting the function of the ATPase by blocking the ion channel (Figure 3C). 
[0052] The model of the binding site allows for two primary observations: First, that the 
consensus binding orientation is highly conserved. With no outliers, the docking modeling 
results in a single binding mode and suggests that only the cardiotonic steroid moiety, lactone 
ring and alpha sugar contribute to binding. Second, while the overall domain of binding is 
very consistent, the individual orientation of the steroid backbones of the inhibitors is not 
uniform indicating that this moiety is acting as a scaffold which provides appropriate spacing 
for the functional groups, rather than as a pharmacophore element essential to binding. 
Furthermore, the number of sugars has no effect on the binding consensus. For example, 
digitoxin with two or three sugars maintains an identical binding orientation (Figure 3D). 
This observation agrees not only with previous in vitro studies describing the orientation of 
the sugars of bound inhibitors, but furthermore corresponds with our own observations based 
on our CoMFA models. 

[0053] One factor that influences the Gibb's free energy (AG) of binding of a small 
molecule to a receptor is the formation of hydrogen bonds between the compound and the 
receptor. One of the most widely used definitions to determine hydrogen bonding potential 
for secondary structure determination is the Database of Secondary Structure in Proteins or 
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DSSP. Systematic measurement of the distance between all potential hydrogen bond donors 
and/or acceptors of the ligands and the residues comprising the binding pocket of the 
invention, the potential residue interactions of each compound within the binding site were 
mapped Of the 16 molecules docked to the la-subunit, 15 had the potential to form 
hydrogen bonds with the side chain or backbone nitrogen of residue D121. Biological data 
has suggested that D121 residue is essential for inhibitor binding as a single point mutation 
(D121N) at this location decreases sensitivity to ouabain by 1000 fold. Four additional 
residues also appear to be involved in molecular inhibition: the side chains of Ql 1 1 and 
N122, and the backbone nitrogens of E908 and M973. The data provided by the methods and 
models of the present invention show that 15 of the 16 molecules docked could potentially 
form hydrogen bonds with QUI, and mutants of the loc-subunit oftheNa, K-ATPase with R 
replacing residue Ql 1 1, or R and D replacing residues Ql 1 1 and N122 respectively, are 8.3 
fold and 1 150 fold less sensitive to ouabain than the wild-type sheep catalytic subunit. In 
addition, all 16 molecules were able to form hydrogen bonds with the backbone nitrogen of 
E908, and 12 molecules were shown to have the potential to interact with M973. 
[0054] The effects of these mutations on the binding of the cardioglycosides can further 
be predicted by modeling the loc-subunit of the Na, K-ATPase incorporating either the 
Ql 1 1R or D121N mutations. In addition, the mutations were introduced in the models of the 
human and rat la-subunits of the Na, K-ATPase whose sequence differs from that of the 
sheep by 23 residues (or 2%) and 35 residues (3.5%) respectively. The human la-subunit 
remains sensitive to ouabain while the rat subunit is insensitive to ouabain inhibition. The 
docking study with ouabain was repeated for both the human and rat la-subunits, and the 
mutant structures maintaining identical parameters as previously. Ouabain binding to the 
human la-subunit of the Na, K-ATPase is identical to that observed in the original docking 
study. The binding affinity of ouabain for both the human and the sheep la-subunit of the 
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Na, K-ATPase is analogous, thus inclusion of the human model acts as a control for the 
following mutational studies. The incorporation of either single point mutation (Ql 1 1R or 
D121N) results in the reorientation of the lactone moiety from H1-H2 toward H3-H4. 
Furthermore, the rat (which contains both Ql 1 1R and N122D point mutations) forces 
ouabain to abandon the consensus orientation observed with the human and sheep subunits. 
[0055] The mutational effects with respect to ouabain sensitivity, described above, can be 
further explored computationally as a decrease in receptor sensitivity to an inhibitor often 
reflects diminished binding affinity. As such, binding energies (AEbmding) associated with 
various ouabain-enzyme interactions resulting from the docking studies were calculated and 
then compared with available inhibition data. In order to quantify the reorientation of ouabain 
with respect to the affinity decreases observed experimentally, the relative change in binding 
energy (AEbmding) were calculated. Values of AE bm dmg were computed as the difference in 
calculated energy between the ouabain-receptor complex (Ecompiex) and the sum of the 
energies calculated separately for the free ligand (Ei igand ) and the free receptor (Er CCCpt0T ). 
[0056] The molecular mechanics calculations included bonded and nonbonded (vdW and 
Coulombic) energy terms and the electrostatic component of solvation. The AE b indin S data 
(Table 3) shows a strong statistical correlation with the experimentally derived inhibitory 
effects of the mutants (Pearson correlation r = 0.81; Spearman correlation r = 0.88). The 
calculated values of AEbmding for the wild-type sheep and human al -subunits are consistent 
with their nearly equivalent ouabain affinity. Furthermore, introduction of the Ql 1 1R 
mutation which is known to lower ouabain's affinity for the sheep enzyme by ~10-fold 
resulted in a corresponding decrease in the calculated binding affinity. Furthermore, the rat 
ccl-subunit, or introduction of the D121N into the sheep al model, yielded AE b i n dmg values 
consistent with the dramatic loss (-1000 fold) loss in ouabain inhibition associated with 
either of these enzymes. 
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[0057] As the results of our molecular modeling studies have implicated two additional 
residues as important components of the digitalis binding site, it was sought to further 
investigate the role of these amino-acid residues with respect to ouabain binding. While in 
vitro mutation studies are beyond the scope of the present study, these mutations were 
simulated in silico to determine the influence of these residues on the affinity of ouabain. 
Mutant sheep al-subunits were constructed incorporating either an E908L or M973A 
mutation and docking analyses with ouabain were performed.. The consensus orientation of 
ouabain selected for analysis was obtained as described above. 

[0058] Interestingly, the consensus orientation observed from the docking ensemble 
differed from the consensus orientation of ouabain bound to the wild-type sheep enzyme 
suggesting that these residues each independently play a role in ouabain binding. Values of 
AEbinding were calculated (Table 4) and compared with the AEbinding values calculated for the 
wild-type and mutant sheep-, rat- and human enzymes. The results indicate that both 
mutants could have a large impact on ouabain binding affinity. 
Table 3. Calculated changes in Binding Energies (AEbinding) 



Species/Mutant* 


Relative Inhibition 


AAEbindiim (kcal/mol) 


Sheep 


1 


-204 


Rat 


1000 


-158 


Human 


1 


-199 


D121N 


1000 


-107 


Q111R 


10 


-167 


E908L 


ukn 


-160 


M973A 


ukn 


-131 



♦Mutation of sheep structure. 



[0059] In another preferred embodiment, the invention provides a method of identifying 
novel ligands of the Na, K-ATPase. Identification of novel ligands can be achieved by 
methods known in the art, such as the screening of libraries and databases for chemicals 
including proteins, peptides, nucleotides, and small molecules that correspond to the 
structural criteria defined by the pharmacophore of the invention. The pharmacophore of the 
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present invention contains the essential hydrophobic and hydrophilic regions for interaction 
(Figure 4). 

[0060] There are essentially three regions of the binding pocket that could be essential for 

binding: a hydrogen bond donor/acceptor corresponding to the potential hydrogen bond 

donor/acceptor of Ql 1 1, an intervening hydrophobic area, and a hydrogen bond donor 

corresponding to the potential hydrogen bond acceptor D121 (represented respectively by 

spheres 1,2 and 3; Figure 4). The two hydrogen bonding groups are 9.95A and 7.63A from 

the hydrophobic core, resulting in a total distance of 17.10 A separating the two hydrogen 

bonding (Figure 4A). In addition to the three pharmacophoric elements discussed above, the 

3D coordinates of E908 (sphere 4; Figure 4) and M973 (sphere 5; Figure 4) provide 

positioning for auxiliary potential hydrogen bond donor/acceptor groups. The rigor of the 

database search criteria can thus be increased in a stepwise fashion by introducing four and 

five point pharmacophores (Table 4). 

Table 4. 
Pharmacophoi 
Three Point 

Four Point 



Five Point 



Sphere 1 


Sphere 2 


Sphere 3 


Sphere 4 


Sphere5 


HBD 


HYD 


HBD 






HBA 


HYD 


HBD 






HBD 


HYD 


HBD 


HBD 




HBA 


HYD 


HBD 


HBA 




HBD 


HYD 


HBD 




HBD 


HBA 


HYD 


HBD 




HBA 


HBD 


HYD 


HBD 


HBA 




HBA 


HYD 


HBD 


HBD 




HBD 


HYD 


HBD 




HBA 


HBA 


HYD 


HBD 




HBD 


HBD 


HYD 


HBD 


HBD 


HBD 


HBD 


HYD 


HBD 


HBD 


HBA 


HBD 


HYD 


HBD 


HBA 


HBD 


HBD 


HYD 


HBD 


HBA 


HBA 


HBA 


HYD 


HBD 


HBD 


HBD 


HBA 


HYD 


HBD 


HBD 


HBA 


HBA 


HYD 


HBD 


HBA 


HBD 


HBA 


HYD 


HBD 


HBA 


HBA 



*HBA represents Hydrogen Bond Acceptor, HBD represents Hydrogen Bond Donor, HYD represents 
Hydrophobic group 
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Distances between points are located in Table 5. 



Table 5. 



Pharmacophore 
Elements 



2 



4 



1 

2 
3 
4 



15 ±0.75 8 ±0.40 18 ±0.90 22 ±1.10 
- 9 ±0.45 6 ±0.30 9 ±0.45 
- 12 ±0.60 14 ±0.70 
7 ±0.35 



* All measurements are in A. 



[0061] However, compounds do not have to satisfy all criteria to be active in inhibiting 
the la subunit of the ATPase. The combination of the initial three pharmacophore elements 
(spheres 1, 2, and 3) collectively represent the minimum 3D criteria necessary to search 
existing databases for novel potential scaffolds. To validate the basic three-point 
pharmacophores of the invention, NCI and Maybridge three-dimensional databases were 
searched for molecules that satisfy the pharmacophore parameters. The pharmacophores of 
the invention allowed the identification of existing cardioglycosides that were independent of 
those used for pharmacophore development in the databases utilizing the pharmacophore. 
Therefore, the three-dimensional selection criteria developed by the present invention mirror 
the physical elements required for inhibitor binding and can be therefore used for the design 
of new therapeutic agents. In addition to the identification of cardioglycosides, the database 
search identified several new potential scaffolds of interest. Based on the models and 
methods of the invention, it was shown that the cardiosteroid backbone of the inhibitors are 
acting merely as hydrophobic spacers, and therefore substitution of scaffolds that maintain 
the three dimensional orientation of the pharmacophore elements should succeed in 
inhibiting the Na, K- ATPase. 

[0062] Thus, in yet another preferred embodiment of the invention, a method of 
designing ligands that will inhibit the Na, K-ATPase is provided. The structural data 
provided by the models and methods of the present invention make it possible to elucidate the 
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mechanism of cardioglycoside inhibition of the la-subunit of the Na, K-ATPase. This in 
turn allows the design of de novo inhibitors of the la-subunit of the Na, K-ATPase. Several 
computer-based schemes using molecular modeling software may be employed to build 
potent and selective inhibitors from these scaffolds. Candidate scaffolds can be placed in the 
binding site of the target molecule, and a ligand can be built by successively bonding other 
fragments to it using methods known in the art. Thus, the Na, K-ATPase model and 
pharmacophores of the present invention will be used to guide the * growth' of a potential 
inhibitor. Public and commercial chemical libraries (ACD, Maybridge, NCI, and WDI) may 
be searched for new molecules or scaffolds that satisfy the query's pharmacophore within 
specified tolerances. Further, computational tools will be employed to match small 
molecules with each other or with receptor topology rapidly based on similarity or 
complementarity in shape and electrostatic properties. 

[0063] Promising candidate inhibitors emerging from the computer-aided design pipeline 
can be evaluated in terms of synthetic feasibility and overall acceptability. Following 
chemical synthesis, binding and inhibition assays known in the art can be used to evaluate the 
in vitro effects of the potential inhibitors. 

Experimental Details 

[0064] The following detailed examples are provided to illustrate the claimed invention 
better and are not to be interpreted as limiting the scope of the invention. To the extent that 
specific materials are mentioned, it is merely for purposes of illustration and is not intended 
to limit the invention. One skilled in the art may develop equivalent means or reactants 
without the exercise of inventive capacity and without departing from the scope of the 
invention. 
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Sequence Alignment between SBRCAla and Na. K-ATPase 

[0065] The sequences of the al-subunits of the human, sheep, and rat Na + , K + -ATPases 
were obtained from GenBank (Accessions nos. P05023, P04074, P06685 respectively). The 
SERCAla crystal structure (1EUL) was obtained from RCSB-PDB and the sequence was 
extracted using Insight II (Accelrys, San Diego, CA). The sequence alignment between 
SERCAla and the al-subunit of the sheep Na + , K + -ATPase was deduced from matrix based 
(PAM250 and Blosum35) approaches fwww.ebi.ac.uk/clustalw/) and by threading the sheep 
sequence to the structural features of SERCAla using THREADER 
(httpr^ioinf.cs.ucl.ac.uk) from University College London and the resulting alignment was 
further manually refined. The amino acid identity resulting from this alignment is 25.6% 
and all local areas of structural and functional importance are well conserved. The 
extracellular gaps resulting from the alignment are localized to loop regions (H1H2, and 
H7H8) and are thus easily accommodated. 

Homology Modeling of the Na. K-ATPase 

[0066] All molecular modeling operations were carried out using Insightll and related 
modules (Accelrys, Inc., San Diego, CA). The three-dimensional coordinates of the 
extracellular and transmembrane domains of the al subunit of sheep Na + JEC + -ATPase were 
assigned. Undefined loops and gaps in the aligned sequences were accommodated by 
random generation of intervening sequences while conserving the integrity of backbone 
orientations in structurally defined, flanking residues. Inclusion of the randomly modeled 
regions required splice site relaxation in order to allow the protein backbone to adopt 
appropriate psi/phi angles. This operation was completed by relaxing the backbone atoms of 
those residues involved in splice sites and subjecting the molecule to a localized energy 
minimization procedure. 
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[0067] After generation of the model, the side chain residues were checked for acceptable 
rotomer positioning using the internal rotomer library and steric clashes were resolved by 
brief energy minimization of the entire molecule (250 iterations, steepest descent). Finally, 
the model was subjected to a full energy minimization through a series of independent steps 
in which hydrogen atoms, side chains, backbone atoms and, ultimately, the entire molecule 
were successively allowed freedom of movement. Similarly, the human and rat sequences 
were aligned to the sheep sequence, and then modeled and refined as described above. Point 
mutations as described herein were incorporated into the fully refined sheep structure, and the 
resulting mutants were again energy minimized through the steps discussed above. This 
procedure implemented the CVFF force field with a distance dependent dielectric function 
(s=Sor, with 8o = 4) until the convergence criterion of 0.04184 kJ/mol (0.01 kcal/mol) change 
in energy between successive iterations was achieved. The dielectric constant € 0 = 4 was 
chosen as it corresponds to that of water in a fast-moving dynamic field. 

Development of the Pharmacophore Model of the Na, K-ATPase 

[0068] The 3D structure of each cardiotonic compound was built using the molecular 
fragment library of Sybyl6.6 (Tripos, Inc., St. Louis, MO) and based on the structure of 
digoxin as determined by X-ray crystallography.(ref) The conformation of the |3- and y-sugars 
was obtained by global energy minimization determinations. As previously described, the 
final structural models were energy minimized using the Tripos force field with the 
Gasteiger-Marsili method for assigning partial atomic charges.(ref) The GOLD operator 
usage probabilities for crossover, mutation, and migration were set to 95%, 95%, and 10%, 
respectively. The maximum distance between hydrogen bond donors and fitting points was 
set to 5A, and the nonbonded cutoff for the van der Waals energies was adjusted to the 
manufacturer recommended settings (10k# = well depth of the van der Waals energy of the 
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atom pair ij). The ligand orientation accepted was the top scoring consensus orientation of 
25 independent genetic algorithm (GA) runs, each with a maximum number of 1000 GA 
operations performed on a single population of 100 individuals. The selection pressure was 
1.1. The centroid of the binding sphere was determined as the intersection of the vectors 
connecting H1H2 with H5H6, H1H2 with H9H10 and H5H6 with H9H10. 

[0069] Binding energy calculations were performed as follows: AEbmding is defined as the 
difference between the potential energy of the ouabain-receptor complex (E CO mpiex) and the 
sum of potential energies of the ouabain (Eijgand) and receptor (Ercceptor). 

AEbinding ~ Ecomplex ~" (Eijgand + Ereceptor) 

Energy values represent the minimum potential energy resulting from the optimization 
procedure described above. Ouabain was allowed full freedom of movement throughout the 
process. Various orientations of ouabain with respect to the ATPase were obtained by 
combining the consensus orientations of ouabain as determined by the GOLD docking, with 
the appropriate homology modeled enzyme. Specifically, the bonded (Stretch, torsion etc.) 
and nonbonded (van der Waals and Columbic terms) energy terms were calculated by 
Discover (InsightlT) whilst the linear Poisson-Boltzmann equation was solved using DelPhi to 
obtain the electrostatic contributions of the solvation energy (reaction field energy). 

Design of Novel Na, K- ATPase Inhibitors 

[0070] Candidate inhibitors identified from the receptor-based pharmacophores provided 
by the present invention will be further screened in terms of bioavailability and toxicological 
profile. Bioavailability is a measure of the value of the hydrophobic/hydrophilic balance, as 
measured by logP, will be predicted using the ClogP program (Daylight Chemical 
Information Systems, Santa Fe, NM). Values of the aqueous solubility and cell membrane 
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permeability will also be predicted using available molecular models in the laboratory as 
additional screens to select only the most promising candidate inhibitors. The toxicological 
profile of potential inhibitors will be determined by calculating relevant endpoints including 
carcinogenicity, Ames mutagenicity, developmental toxicity potential, skin irritancy, ocular 
irritancy, and aerobic biodegradability will be predicted for the candidate inhibitors using the 
QSAR-based TOPKAT program (Accelrys, Inc., San Diego, CA). 

[0071] Several computer-based schemes may be employed to build potent and selective 
inhibitors from these scaffolds. First, a method known as "molecule growth" may be 
employed to place a 'seed' fragment or scaffold in the binding site of the target molecule, and 
to build a ligand by successively bonding other fragments to it using methods such as Small 
Molecule Growth (SmoG) (DeWhitte and Shakhnovich, J. Am. Chem. Soc. 118: 11733- 
11744, 1996)) or GrowMol (Bohacek andMcMartin, J. Am. Chem. Soc. 116: 5560-5571, 
1994). The Na, K-ATPase homology model of the present invention will be used to guide the 
'growth' of a potential inhibitor. 

[0072] Second, public and commercial chemical libraries (ACD, Maybridge, NCI, and 
WDI) may be searched using a fragment-based development method commonly known as 
"sub-structure searching." This method uses a 3D pharmacophore as a query to search a 
chemical library for new molecules (or scaffolds) that satisfy the query's pharmacophore 
within specified tolerances. Third, a computational tool will be employed as a ray-tracing 
scheme to rapidly match small molecules with each other (or even with receptor pockets) 
based on similarity (or complementarity) in shape and electrostatic properties. It is a 
powerful and efficient method for drug discovery and particularly for the discovery of new 
scaffolds that share the same properties but not necessarily the same structural formula as the 
query molecule. 
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[0073] The design strategies used for chemical synthesis will depend on the structures of 
the proposed inhibitors. Promising candidate inhibitors emerging from the computer-aided 
design pipeline will be evaluated in terms of synthetic feasibility and overall acceptability. 
Following chemical synthesis, binding and inhibition assays will be used to evaluate the in 
vitro effects of the potential inhibitors. Assays will be conducted as describe previously 
(Biochemistry 2002, 41, 1 137-1 148). Briefly Na, K-ATPase will be purified from the outer 
medulla of frozen lamb kidney. Competitive binding assays will be performed in the 
presence and absence of varying concentrations of the potential inhibitors. Inhibition studies 
will be carried out with and without varying concentrations of the potential inhibitors, and the 
NADPH oxidation rate in the presence and absence of the inhibitors will be determined. 
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CLAIMS 

What is claimed is: 

1 . A novel pharmacophore model as defined by the parameters of Table 4 and Table 

5. 

2. The novel pharmacophore model of claim 1, wherein scaffold molecules derived 
therefrom can be used as a basis for compounds directed to inotropic Na, K-ATPase 
inhibition. 

3. The novel pharmacophore model of claim 1, wherein the model produces an Na, K- 
ATPase inhibitor compound of the formula: 



o 




wherein Rl, R2, R3 and R4 can be any organic functional group containing a hydrogen bond 
donor or a hydrogen bond acceptor and X is any element or group that allows the compound 
to retain inotropic activity. 

4. The novel pharmacophore model of claim 3, wherein X is N, O, S, or C. 

5. A method of using a pharmacophore model to create an Na, K-ATPase 
inhibitory compound comprising the steps of: 
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(a) creating alignment between SERCA andNa, K-ATPase, wherein SERCA is a template; 

(b) transferring coordinates from the template to a model for structurally conserved regions; 

(c) generating variable regions; 

(d) refining the model through energy minimization steps; and 

(e) performing docking analysis of prospective drug candidates. 

6. The method of claim 5, further comprising the steps of: 

(f) delineating the essential phannacophoric elements for high binding affinity; 

(g) searching databases of known compounds using the restraints as implicated by the 
pharmacophore with allowable tolerances; and 

(h) utilizing de novo rational drug design and computer aided molecular modeling to design 
novel compounds using the restraints as implicated by the pharmacophore with allowable 
tolerances. 

7. The method of claim 6, wherein the allowable tolerances in steps (g) and (h) is 

±10%. 

8. The method of claim 5, wherein step (a) is comprised of dynamic 
programming and threading. 

9. The method of claim 5, wherein SERCA is SERCA1 a. 

10. The method of claim 5, wherein the steps are carried out using a computer- 
readable medium having computer-executable instructions. 

1 1 . The method of claim 1 0, wherein the steps are carried out using molecular 
modeling software. 

12. A method of treating an individual with a heart disease comprising 
administering a therapeutically effective amount of an novel inotropic compound created 
using a novel pharmacophore model as defined by Table 4 and Table 5. 
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13. The method of claim 12, wherein the novel pharmacophore model produces 
novel inotropic drugs of the formula of: 




14. The method of claim 13, wherein the novel drugs have a wider therapeutic index than 
either ouabain or digoxin. 

15. The method of claim 12, wherein the heart disease treated is congestive heart failure and 
supraventricular arrhythmia. 

16. The method of claim 12, wherein the novel inotropic compound is administered in a 
pharmaceutical^ acceptable carrier. 

17. The method of claim 12, wherein the novel inotropic compound is administered 
parenterally or orally. 

18. The method of claim 12, wherein residues Ql 1 1, D121, E908 and M973 are unaltered. 
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Figure 3B 
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SEQUENCE LISTING 
<110> University of Medicine and Dentistry of New Jersey 

<120> NOVEL PHARMACOPHORE FOR THE DISCOVERY AND TESTING OF NA, K-ATPASE 
INHIBITOR COMPOSITIONS AND METHODS FOR THEIR USE IN TREATING 
CARDIOVASCULAR DISEASES AND CONDITIONS 

<130> 54704. 8061. WO00 

<150> 60/425,037 
<151> 2002-11-07 

<160> 2 

<170> Patentln version 3.2 

<210> 1 

<211> 1001 

<212> PRT 

<213> Homo sapiens 

<400> 1 

Met Glu Ala Ala His Ser Lys Ser Thr Glu Glu Cys Leu Ala Tyr Phe 
15 10 15 



Gly Val Ser Glu Thr Thr Gly Leu Thr Pro Asp Gin Val Lys Arg His 
20 25 30 



Leu Glu Lys Tyr Gly His Asn Glu Leu Pro Ala Glu Glu Gly Lys Ser 
35 40 45 



Leu Trp Glu Leu Val lie Glu Gin Phe Glu Asp Leu Leu Val Arg lie 
50 55 ~ 60 



Leu Leu Leu Ala Ala Cys lie Ser Phe Val Leu Ala Trp Phe Glu Glu 
65 70 75 * 80 



Gly Glu Glu Thr lie Thr Ala Phe Val Glu Pro Phe Val He Leu Leu 
85 90 95 



He Leu He Ala Asn Ala He Val Gly Val Trp Gin Glu Arg Asn Ala 
100 105 110 



Glu Asn Ala He Glu Ala Leu Lys Glu Tyr Glu Pro Glu Met Gly Lys 
115 120 125 



Val Tyr Arg Ala Asp Arg Lys Ser Val Gin Arg He Lys Ala Arg Asp 
130 135 140 
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lie Val Pro Gly Asp lie Val Glu Val Ala Val Gly Asp Lys Val Pro 
145 150 155 160 



Ala Asp lie Arg lie Leu Ser lie Lys Ser Thr Thr Leu Arg Val Asp 
165 170 175 



Gin Ser He Leu Thr Gly Glu Ser Val Ser Val He Lys His Thr Glu 
180 185 190 



Pro Val Pro Asp Pro Arg Ala Val Asn Gin Asp Lys Lys Asn Met Leu 
195 200 205 



Phe Ser Gly Thr Asn He Ala Ala Gly Lys Ala Leu Gly He Val Ala 
210 215 220 



Thr Thr Gly Val Ser Thr Glu He Gly Lys He Arg Asp Gin Met Ala 
225 230 235 240 



Ala Thr Glu Gin Asp Lys Thr Pro Leu Gin Gin Lys Leu Asp Glu Phe 
245 250 255 



Gly Glu Gin Leu Ser Lys Val He Ser Leu He Cys Val Ala Val Trp 
260 265 270 



Leu He Asn He Gly His Phe Asn Asp Pro Val His Gly Gly Ser Trp 
275 280 285 



He Arg Gly Ala He Tyr Tyr Phe Lys He Ala Val Ala Leu Ala Val 
290 295 300 



Ala Ala He Pro Glu Gly Leu Pro Ala Val He Thr Thr Cys Leu Ala 
305 310 315 320 



Leu Gly Thr Arg Arg Met Ala Lys Lys Asn Ala He Val Arg Ser Leu 
325 330 335 



Pro Ser Val Glu Thr Leu Gly Cys Thr Ser Val He Cys Ser Asp Lys 
340 345 350 



Thr Gly Thr Leu Thr Thr Asn Gin Met Ser Val Cys Lys Met Phe He 
355 360 365 



He Asp Lys Val Asp Gly Asp Phe Cys Ser Leu Asn Glu Phe Ser He 
370 375 380 
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Thr Gly Ser Thr Tyr Ala Pro Glu Gly Glu Val Leu Lys Asn Asp Lys 
385 390 395 400 



Pro lie Arg Ser Gly Gin Phe Asp Gly Leu Val Glu Leu Ala Thr lie 
405 410 415 



Cys Ala Leu Cys Asn Asp Ser Ser Leu Asp Phe Asn Glu Thr Lys Gly 
420 425 430 



Val Tyr Glu Lys Val Gly Glu Ala Thr Glu Thr Ala Leu Thr Thr Leu 
435 440 445 



Val Glu Lys Met Asn Val Phe Asn Thr Glu Val Arg Asn Leu Ser Lys 
450 455 460 



Val Glu Arg Ala Asn Ala Cys Asn Ser Val lie Arg Gin Leu Met Lys 
465 470 475 480 



Lys Glu Phe Thr Leu Glu Phe Ser Arg Asp Arg Lys Ser Met Ser Val 
485 490 495 



Tyr Cys Ser Pro Ala Lys Ser Ser Arg Ala Ala Val Gly Asn Lys Met 
500 505 510 



Phe Val Lys Gly Ala Pro Glu Gly Val He Asp Arg Cys Asn Tyr Val 
515 520 525 



Arg Val Gly Thr Thr Arg Val Pro Met Thr Gly Pro Val Lys Glu Lys 
530 535 540 



He Leu Ser Val He Lys Glu Trp Gly Thr Gly Arg Asp Thr Leu Arg 
545 550 555 560 



Cys Leu Ala Leu Ala Thr Arg Asp Thr Pro Pro Lys Arg Glu Glu Met 
565 570 575 



Val Leu Asp Asp Ser Ser Arg Phe Met Glu Tyr Glu Thr Asp Leu Thr 
580 585 590 



Phe Val Gly Val Val Gly Met Leu Asp Pro Pro Arg Lys Glu Val Met 
595 600 605 



Gly Ser He Gin Leu Cys Arg Asp Ala Gly He Arg Val He Met He 
610 615 620 
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Thr Gly Asp Asn Lys Gly Thr Ala lie Ala lie Cys Arg Arg He Gly 
625 630 635 640 



He Phe Gly Glu Asn Glu Glu Val Ala Asp Arg Ala Tyr Thr Gly Arg 
645 650 655 



Glu Phe Asp Asp Leu Pro Leu Ala Glu Gin Arg Glu Ala Cys Arg Arg 
660 665 670 



Ala Cys Cys Phe Ala Arg Val Glu Pro Ser His Lys Ser Lys He Val 
675 680 685 



Glu Tyr Leu Gin Ser Tyr Asp Glu He Thr Ala Met Thr Gly Asp Gly 
690 695 700 



Val Asn Asp Ala Pro Ala Leu Lys Lys Ala Glu He Gly He Ala Met 
705 710 715 720 



Gly Ser Gly Thr Ala Val Ala Lys Thr Ala Ser Glu Met Val Leu Ala 
725 730 735 



Asp Asp Asn Phe Ser Thr He Val Ala Ala Val Glu Glu Gly Arg Ala 
740 745 750 



He Tyr Asn Asn Met Lys Gin Phe He Arg Tyr Leu He Ser Ser Asn 
755 760 765 



Val Gly Glu Val Val Cys He Phe Leu Thr Ala Ala Leu Gly Leu Pro 
770 775 780 



Glu Ala Leu He Pro Val Gin Leu Leu Trp Val Asn Leu Val Thr Asp 
785 790 795 800 



Gly Leu Pro Ala Thr Ala Leu Gly Phe Asn Pro Pro Asp Leu Asp He 
805 810 ~ 815 



Met Asp Arg Pro Pro Arg Ser Pro Lys Glu Pro Leu He Ser Gly Trp 
820 825 830 



Leu Phe Phe Arg Tyr Met Ala He Gly Gly Tyr Val Gly Ala Ala Thr 
835 840 845 



Val Gly Ala Ala Ala Trp Trp Phe Met Tyr Ala Glu Asp Gly Pro Gly 
850 855 ~ 860 
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Val Thr Tyr Hi s 1 Gin Leu Thr His Phe Met Gin Cys Thr Glu Asp His 
865 870 875 880 



Pro His Phe Glu Gly Leu Asp Cys Glu He Phe Glu Ala Pro Glu Pro 
885 890 895 



Met Thr Met Ala Leu Ser Val Leu Val Thr He Glu Met Cys Asn Ala 
900 905 910 



Leu Asn Ser Leu Ser Glu Asn Gin Ser Leu Met Arg Met Pro Pro Trp 
915 920 925 



Val Asn He Trp Leu Leu Gly Ser He Cys Leu Ser Met Ser Leu His 
930 935 940 



Phe Leu lie Leu Tyr Val Asp Pro Leu Pro Met He Phe Lys Leu Lys 
945 950 955 960 



Ala Leu Asp Leu Thr Gin Trp Leu Met Val Leu Lys He Ser Leu Pro 
965 970 975 



Val He Gly Leu Asp Glu He Leu Lys Phe He Ala Arg Asn Tyr Leu 
980 985 990 



Glu Asp Pro Glu Asp Glu Arg Arg Lys 
995 * 1000 



<210> 2 

<211> 1016 

<212> PRT 

<213> Ovis sp. 

<400> 2 

Gly Arg Asp Lys Tyr Glu Pro Ala Ala Val Ser Glu His Gly Asp Lys 
15 10 15 



Lys Lys Ala Lys Lys Glu Arg Asp Met Asp Glu Leu Lys Lys Glu Val 
20 " 25 30 



Ser Met Asp Asp His Lys Leu Ser Leu Asp Glu Leu His Arg Lys Tyr 
35 40 45 



Gly Thr ABp Leu Asn Arg Gly Leu Thr Thr Ala Arg Ala Ala Glu He 
50 55 60 
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Leu Ala Arg Asp Gly Pro Asn Ala Leu Thr Pro Pro Pro Thr Thr Pro 
65 70 75 80 



Glu Trp Val Lys Phe Cys Arg Gin Leu Phe Gly Gly Phe Ser Met Leu 
85 90 95 



Leu Trp lie Gly Ala Val Leu Cys Phe Leu Ala Tyr Gly He Gin Ala 
100 105 110 



Ala Thr Glu Glu Glu Pro Gin Asn Asp Asn Leu Tyr Leu Gly Val Val 
115 120 125 



Leu Ser Ala Val Val He He Thr Gly Cys Phe Ser Tyr Tyr Gin Glu 
130 135 140 



Ala Lys Ser Ser Lys He Met Glu Ser Phe Lys Asn Met Val Pro Gin 
145 150 155 160 



Gin Ala Leu Val He Arg Asn Gly Glu Lys Met Ser He Asn Ala Glu 
165 170 175 



Glu Val Val Val Gly Asp Leu Val Glu Val Lys Gly Gly Asp Arg He 
180 185 190 



Pro Ala Asp Leu Arg lie lie Ser Ala Asn Gly Cys Lys Val Asp Asn 
195 200 205 



Ser Ser Leu Thr Gly Glu Ser Glu Pro Gin Thr Arg Ser Pro Asp Phe 
210 215 220 



Thr Asn Glu Asn Pro Leu Glu Thr Arg Asn He Ala Phe Phe Ser Thr 
225 230 235 240 



Asn Cys Val Glu Gly Thr Ala Arg Gly lie Val Val Tyr Thr Gly Asp 
245 250 255 



Arg Thr Val Met Gly Arg lie Ala Thr Leu Ala Ser Gly Leu Glu Gly 
260 265 270 



Gly Gin Thr Pro He Ala Ala Glu He Glu His Phe He His He He 
275 280 285 



Thr Gly Val Ala Val Phe Leu Gly Val Ser Phe Phe lie Leu Ser Leu 
290 295 300 
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lie Leu Glu Tyr Thr Trp Leu Glu Ala Val lie Phe Leu lie Gly lie 
305 310 315 320 



He Val Ala Asn Val Pro Glu Gly Leu Leu Ala Thr Val Thr Val Cys 
325 330 335 



Leu Thr Leu Thr Ala Lys Arg Met Ala Arg Lys Asn Cys Leu Val Lys 
340 345 350 



Asn Leu Glu Ala Val Glu Thr Leu Gly Ser Thr Ser Thr He Cys Ser 
355 360 365 



Asp Lys Thr Gly Thr Leu Thr Gin Asn Arg Met Thr Val Ala His Met 
370 375 380 



Trp Phe Asp Asn Gin He His Glu Ala Asp Thr Thr Glu Asn Gin Ser 
385 390 395 400 



Gly Val Ser Phe Asp Lys Thr Ser Ala Thr Trp Leu Ala Leu Ser Arg 
405 410 415 



He Ala Gly Leu Cys Asn Arg Ala Val Phe Gin Ala Asn Gin Asp Asn 
420 425 430 



Leu Pro He Leu Lys Arg Ala Val Ala Gly Asp Ala Ser Glu Ser Ala 
435 440 445 



Leu Leu Lys Cys He Glu Val Cys Cys Gly Ser Val Lys Glu Met Arg 
450 " 455 460 



Glu Arg Tyr Ala Lys He Val Glu He Pro Phe Asn Ser Thr Asn Lys 
465 470 475 480 



Tyr Gin Leu Ser He His Lys Asn Ala Asn Ala Gly Glu Pro Arg His 
485 490 495 



Leu Leu Val Met Lys Gly Ala Pro Glu Arg He Leu Asp Arg Cys Ser 
500 505 510 



Ser He Leu He His Gly Lys Glu Gin Pro Leu Asp Glu Glu Leu Lys 
515 520 525 



Asp Ala Phe Gin Asn Ala Tyr Leu Glu Leu Gly Gly Leu Gly Glu Arg 
530 535 540 
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Val Leu Gly Phe Cys His Leu Met Leu Pro Asp Glu Gin Phe Pro Glu 
545 550 555 560 



Gly Phe Gin Phe Asp Thr Asp Asp Val Asn Phe Pro Val Asp Asn Leu 
565 570 575 



Cys Phe Val Gly Leu lie Ser Met lie Asp Pro Pro Arg Ala Ala Val 
580 585 " 590 



Pro Asp Ala Val Gly Lys Cys Arg Ser Ala Gly He Lys Val He Met 
595 600 605 



Val Thr Gly Asp His Pro He Thr Ala Lys Ala He Ala Lys Gly Val 
610 615 620 



Gly He He Ser Glu Gly Asn Glu Thr Val Glu Asp He Ala Ala Arg 
625 630 635 640 



Leu Asn He Pro Val Ser Gin Val Asn Pro Arg Asp Ala Arg Ala Cys 
645 650 655 



Val Val His Gly Ser Asp Leu Lys Asp Met Thr Pro Glu Gin Leu Asp 
660 665 670 



Asp He Leu Lys Tyr His Thr Glu He Val Phe Ala Arg Thr Ser Pro 
675 680 685 



Gin Gin Lys Leu He He Val Glu Gly Cys Gin Arg Gin Gly Ala He 
690 695 700 



Val Ala Val Thr Gly Asp Gly Val Asn Asp Ser Pro Ala Leu Lys Lys 
705 710 715 720 



Ala Asp He Gly Val Ala Met Gly He Ala Gly Ser Asp Val Ser Lys 
725 730 735 



Gin Ala Ala Asp Met He Leu Leu Asp Asp Asn Phe Ala Ser He Val 
740 745 750 



Thr Gly Val Glu Glu Gly Arg Leu He Phe Asp Asn Leu Lys Lys Ser 
755 760 765 



He Ala Tyr Thr Leu Thr Ser Asn He Pro Glu He Thr Pro Phe Leu 
770 775 780 
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lie Phe lie He Ala Asn He Pro Leu Pro Leu Gly Thr Val Thr He 
785 790 795 800 



Leu Cys lie Asp Leu Gly Thr Asp Met Val Pro Ala He Ser Leu Ala 
805 810 815 



Tyr Glu Gin Ala Glu Ser Asp He Met Lys Arg Gin Pro Arg Asn Pro 
820 825 830 



Gin Thr Asp Lys Leu Val Asn Glu Arg Leu He Ser Met Ala Tyr Gly 
835 840 845 



Gin He Gly Met He Gin Ala Leu Gly Gly Phe Phe Thr Tyr Phe Val 
850 855 860 



He Met Ala Glu Asn Gly Phe Leu Pro Asn His Leu Leu Gly He Arg 
865 870 875 880 



Val Thr Trp Asp Asp Arg Trp He Asn Asp Val Glu Asp Ser Tyr Gly 
885 890 895 



Gin Gin Trp Thr Tyr Glu Gin Arg Lys He Val Glu Phe Thr Cys His 
900 905 910 



Thr Ala Phe Phe Val Ser He Val Val Val Gin Trp Ala Asp Leu Val 
915 920 925 



He Cys Lys Thr Arg Arg Asn Ser Val Phe Gin Gin Gly Met Lys Asn 
930 935 940 



Lys He Leu He Phe Gly Leu Phe Glu Glu Thr Ala Leu Ala Ala Phe 
945 950 955 960 



Leu Ser Tyr Cys Pro Gly Met Gly Val Ala Leu Arg Met Tyr Pro Leu 
965 970 "* 975 



Lys Pro Thr Trp Trp Phe Cys Ala Phe Pro Tyr Ser Leu Leu He Phe 
980 985 990 



Val Tyr Asp Glu Val Arg Lys Leu He He Arg Arg Arg Pro Gly Gly 
995 1000 1005 



Trp Val Glu Lys Glu Thr Tyr Tyr 
1010 1015 
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